Oscillation in microRNA Feedback Loop 
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The dynamic behaviors of microRNA and mRNA under external stress are studied with biological 
experiments and mathematics models. In this study, we developed a mathematic model to describe 
the biological phenomenon and for the first time reported that, as responses to external stress, the 
expression levels of microRNA and mRNA sustained oscillation. And the period of the oscillation 
is much shorter than several reported transcriptional regulation negative feedback loop. 
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Since the first microRNA (miRNA) lin-4 was discov- 
ered two decades ago, microRNA has been a rising star in 
the field of biology. It has been proved to be related with 
lots of essential processes in cells, such as DNA damage 
response, cell-cycle checkpoints, programmed cell death 
and metabolism [2[. As known to all, microRNA executes 
its functions as a regulator by binding to complementary 
sequences on target messenger RNA transcripts (mRNA) 
and forming RNA-induce silencing complex (RISC), after 
transcription but before translation. Length of a mature 
miRNA is just about 22 nt and a gene-coding mRNA is 
much bigger and needs to translate into a protein. There- 
fore, comparing with the regulation by a protein, the reg- 
ulatory function of a miRNA is a more economical means 

In cell signaling pathway, negative feedback loop is a 
very normal motif, consisted of not only proteins but 
also other non-coding RNAs including miRNAs. The 
advantage of this topological mode lies in its stability 
|3| . There are many researches on protein-protein neg- 
ative feedback loop, e.g., mdm2 and p53 [J-|6|, but few 
on miRNA regulations |7H1Q|. Based on the special bi- 
ological functions of miRNAs, we built a stress-induced 
miRNA-protein dynamical model which contains a neg- 
ative feedback loop. Then we designed a series of exper- 
iments to observe dynamic behaviors of this regulation 
on molecular level. We made a first report on miRNA 
oscillation in cells exposed to radiation and these results 
generally support our model. Furthermore, this model 
could be used to construct dynamical networks of pro- 
teins and other non-coding RNAs. 

A scheme of miRNA feedback loop is shown in Fig. 
[T^a). The mRNA transcription of a certain gene increases 
under stress, and its corresponding protein, as a tran- 
scriptional factor, actives the expression of the miRNA 
which represses the transcription of the mRNA in return. 
A three- variable ODEs is considered to describe the sys- 
tem, 
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FIG. 1: (a) A scheme of miRNA feedback loop. A signal 
of radiation induces transcription of mRNA, and its corre- 
sponding protein, as a transcriptional factor, active the tran- 
scription miRNA. Then the miRNA inhibits the expression of 
protein by inducing the degradation of mRNA. (b) A scheme 
of p53-mdm2 negative feedback loop. 



where [mRNA], [protein], [miRNA] are the expression 
levels of the mRNA, protein, and miRNA, respectively. 
And Sig is the signal that stimulates the responding ef- 
fects of cells and activates the transcription of mRNA. 
For simplicity, Sig is just a constant here. As will be 
shown below, a non-zero Signal is a necessary condition of 
the excitation of a sustained oscillation. d m RNA, d prot eim 
and d m iRNA are degradation rates of the three molecules. 
kmRNA-protein and kp rot ein-miRNA are the rates of pro- 
tein translation and miRNA transcription. Only a linear 
model is considered. Finally, for the regulation of miRNA 
and mRNA, as RISC is a stable machine [2j, we tend to 
consider the condition that miRNA could be recycled, 
and then the Michael-Mathen function is considered. By 
renormalizing time and parameters according to 
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The diurnal oscillation of miRNA and mRNA are re- 
ported in recent works [11]. In this paper, we report that 
this phenomenon can also be observed in stress-response 
process. And the dynamical behaviors of Eq. (J3j) with 
proper parameters shows sustained oscillation, which can 
be used to describe the biological phenomenon. Fig. [2] 
shows an example. Moreover, it can also be observed that 
the periods of the three oscillation curves are almost the 
same, but their phases are different. The oscillation of 
mRNA is firstly excited by a non-zero signal, and fol- 
lowed by the oscillation of the protein. And then the 
expression level of miRNA also rises and begins to re- 
press the production of mRNA. Thus the oscillations of 
mRNA and miRNA are almost anti-phase. 

Since there have been many studies on transcription, 
translation and the degradation of mRNA and protein, 
and miRNA regulation negative feedback loop can be 
analogous to the well studied p53-mdm2 feedback loop 
shown in Fig. DJb), the parameters about these pro- 
cesses are easy to be chosen [12j. But the transcription 
and degradation rates of miRNA and the parameters of 
miRNA- mRNA interaction are still unclear. Thus, these 
parameters will be mainly focused on in this paper. Since 
the degradation rate of miRNA (d m i) has been reduced, 
only are three parameters (k mirri , k pmi , j m im) left. 

Because all parameters are positive, system (J3]) has 
only one positive fixed point. And as will be shown below, 
~ jmim < < kmim , kpmi • Thus the fixed point can be 
written by 
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To assess whether Eq. (J3]) can display oscillations, we 
have searched for parameter values where the fixed point 
specified by @ loses stability to a periodic solution via 
a Hopf bifurcation, that is, a pair of conjugate eigenval- 
ues of the linearized problem cross the imaginary axis. 
Because Eq. ((3]) has three variables, the characteristic 
polynomial of Jacobi matrix in the fixed point (j4]) is a 
cubic monic polynomial, which is denoted by 
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FIG. 2: The time course of the expression of mRNA (green 
line), protein (red line) and miRNA (blue line) described by 
Eq. © with S = 0.9, dm = 0.8, kmim = 6, jmim = 0.001, 
dp — 0.8, kmp — 0.8 dmi — 3, kmpi — 10, and mi — 0.135, 
m = 0.047, p = 0.028 at t = 0. 



The roots of Eq. ([5]) are the eigenvalues of the linearized 
form of Eq. ([3]). They describe the dynamics behavior 
of Eq. (|3|) around the fixed point. If the determinant of 
Eq. © 
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is positive, the linearized Eq. (|3]) will has one real 
eigenvalue (denoted by Ai) and two conjugate imaginary 
eigenvalues (A2 and A3). That is, the oscillation of system 
Eq. (J3j) will be excited around the fixed point Eq. ((4]). 
But, a sustained oscillation requires that the real parts 
of A2 and A3 cross 0, i.e., the Hopf bifurcation. Thus 
according to the general formula of the roots of Eq. (J5]), 
the critical point of Hopf bifurcation can be written as 
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where T = bc/6 - b 3 /27 - d/2 and K = c/3- b 2 /9. Then, 
to cubing both sides of Eq. ©, it will be 
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Finally, we have 



R = d-bc = 0. 



(8) 



(9) 



That is, the condition of the existence of a sustained os- 
cillation in Eq. (J3j) is 



R>0. 
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Thus, to substitute the parameters of Eq. (J3]) into Eq. 
(|6j) and Eq. (fTD|) and to solve the equation R = 0,one can 
obtain the phase graph in parameter space (kmim, jmim, 

fcpmi ) ' 

It is easy to be tested that, no matter in k m im-kpmi 
plane or in k m im-jmim plane, the determinant Eq. (|6j) has 




FIG. 3: The attenuate oscillation area (labeled by A) and the 
sustained oscillation area (labeled by S) in (a) k m im-k p mi 
and (b) jmim-k m im parameter space, which are divided by 
the critical line satisfied R = 0. 





FIG. 5: The experimental data of Dicer mRNA (green line), 
Dicer protein (red line) and miRNA-3928 (blue line) and the 
analysis results of model Eq. Q. 



FIG. 4: (a) The rm curves of Eq. © with 1000 different 
parameter values, (b) rrii, ra, p curves of Eq. J3} with S = 0.9, 
dm = 0.8,kmim = 6, jmim = 0.001, dp = 0.8, kmp = 0.8701 
dmi — 3, kmpi — 7.3298, and mi — 0.135, m — 0.047, p — 
0.028 at t = 0. (c) The curves of < m >, < p > and < rm > 



positive values in a broad parameter area. That means 
an oscillation solution of the system (|3J) is usual and easy 
to be obtained. But, a sustained oscillation is unusual. 
Firstly, the oscillation parameter area(i? > 0) in k m im- 
kpmi plane is shown in Fig. [3(a). This picture can be ex- 
plained in the sense that an over faster transcription rate 
of miRNA (k pm i is big) and an over powerful repression 
to mRNA (kmim is big simultaneously) will dramatically 
reduce the production of mRNA and protein, and then 
oscillation attenuates. Furthermore, in k m im-jmim plane 
(Fig. |3fb)), R > can only be achieved when jmim is 
very small. That is, Eq. ([3]) is a weak nonlinear dynam- 
ics system. Since the Michael-Mathen coefficient jmim is 
the ratio of the rate of a reaction and that of its inverse 
reaction, a very small jmim means the reaction is basi- 
cally irreversible. And this feature can be observed in 
miRNA-mRNA interaction procession. 

Empirically speaking, the oscillation of miRNA and 
protein in a single cell will not reduce under a constant 
drive, but the peak could be broadened because of the 
non-synchronization among the cells. The model Eq. ((3]) 
only describes mono cell behavior. But in real biolog- 
ical systems, the rates of transcription, translation and 
other regulations vary in different cells. This will lead to 
different parameters when the model Eq. ((3]) is used to 
describe different cells. Thus, when the problem refers 



to the behaviors of millions of cells, the parameter ro- 
bustness of the system Eq. (|3j) must be discussed. Fig. 
[3(a) shows the curves of the system Eq. (|3|) with differ- 
ent parameters. It can be found that the system Eq. (J3J) 
is parameter robust. Even though the environments can 
affect some specific biological processes, the existence of 
sustained oscillation will not be dramatically changed by 
small perturbation of environments. On the other hand, 
the amplitudes and periods of oscillations depend on pa- 
rameters values. Therefore, the collective behavior of 
cells will show amplitude attenuation. To illustrate this 
picture, Eq. ([3]) is simulated by using N= 10000 differ- 
ent parameter values (the expectation of these parameter 
values is equal to that we used before and their standard 
deviation is less than 30% of the values that we used 
before). Then the averages over the 10000 dynamical 
processes are calculated: < m(t) >= ^2m k (t)/ 10000, 

< P(fi) >= $> fc (t)/10000, < rrn(t) >= £raf (t)/10000, 
where k is the index of different parameter values. As 
shown in Fig. E(b) an d Fig. E(c), the oscillation of 

< m >, < p > and < rrii > attenuate. 

To testify the theoretical result, an experiment is de- 
signed. We chose human cervical carcinoma HeLa cell 
line for the biological experiment. Dicer protein is re- 
quired to process miRNA precursor to mature miRNAs, 
and meanwhile is a target of miRNA-3928 . So Dicer and 
miRNA-3928 form a negative feedback loop as shown in 
Fig. [1(a). HeLa cells were cultured in RPMI-1640 media 
and exposed to 2 Gy X-rays at room temperature. The 
qRT-PCR (quantitative real-time reverse transcription- 
PCR) is used to measure time-sequencing expression lev- 
els of miRNA-3928 and Dicer mRNA. After irradiation, 
DNA repair were activated in cells. Unless the damages 
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FIG. 6: Western blotting picture of Dicer protein till 2.25 
hours after radiation, in which the oscillation of its concen- 
tration can be found. And the oscillation period is about 1 
hour. 



are repaired, cells would not reenter cell cycle. Under the 
conditions of dose used in this experiment, the duration 
of repair is about 6-8 hours. The drive would disappear at 
that time [14[ . The results are shown by the circle points 
in Fig. [5] And Western blotting is used to measure the 
expression level of Dicer protein. The concentration of 
Dicer is shown by circle points in Fig. [5j and the oscilla- 
tion of Dicer protein can also be found in one of Western 
blotting pictures shown in Fig. [6j To get the numeri- 
cal data, Gray value analysis is used to process pictures. 
And it can be found that the concentrations of miRNA- 
3928, Dicer and its mRNA show an obvious oscillation 
behavior after radiation. The oscillation period is about 
1 hour, and the amplitudes of the oscillations attenu- 
ate after several oscillation periods. The phases of these 
oscillations are different, opposite change trend between 
the concentration of Dicer and that of miRNA-3928 can 
always be observed. 



On the other hand, the transcriptional regulation is 
also popular in biological processes. The oscillations of 
proteins p53 and mdm2 with a phase differences can also 
be observed and various mathematical models were con- 
structed to explain the experiments. But a significant 
difference between transcriptional regulation and miRNA 
regulation is the oscillation periods. The period of p53 
and mdm2 oscillation is about 6 hours, while that of 
miRNA-3928 and Dicer oscillation is only 1 hour. This 
shows that the response time of miRNA regulation is far 
more fasten than that of transcriptional factor regula- 
tion. Obviously, a faster response is very important to 
radiation response, particularly to the repairs of radia- 
tion damages. 

In conclusion, we build a model to describe the 
miRNA-protein negative feedback loop, which is tested 
by the biology experiment of Dicer and miRNA-3928 reg- 
ulation loop. Through the numerical simulation and ex- 
periments, the oscillations of the expression of mRNA, 
protein and miRNA and their phase difference are ob- 
served. We find that the response time of miRNA reg- 
ulation is far faster than that of the transcriptional fac- 
tor regulation. The results of simulations are consistent 
with those of experiments well. These results provide 
estimates of the effective biochemical parameter in the 
processes involving miRNA and could be used to con- 
struct dynamical networks of proteins and other non- 
coding RNAs. 
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